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Abstract 

Supplementary information for our manuscript, entitled 'Spontaneous Skyrmion Ground States 



of Magnetic Metals', cond-mat/0603103 , is presented. The physical nature of the gradient terms 
of our generalized micromagnetic model for ferromagnets with softened longitudinal fluctuations 
is explained. The relationship of our micromagnetic model with the spin fluctuation theory of 
itinerant-electron magnets is discussed. Experimental estimates of the parameter n, which accounts 
for an effective reduced longitudinal stiffness, are presented for real materials from published polar- 
ized neutron scattering experiments on EuS, Ni and MnSi. The available experimental data clearly 
show that n is significantly reduced for the latter two systems. It is suggested that particle-hole 
excitations are at the root of this longitudinal softness in itinerant-electron ferromagnets. The cur- 
rent status of the experimental evidence supporting spontaneous, amorphous skyrmion textures in 
MnSi and other materials is reviewed. Finally, we also address the general potential of skyrmion 
textures in chiral magnets for other fields of physics. 
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I. OUTLINE 



An important ingredient for the prediction of spontaneous skyrmion ground states in 



magnetic metals with chiral interactions in our manuscript cond-mat/0603103 |l| is the 
assumption that the longitudinal magnetic stiffness is particularly soft. In section |H] we 
consider the theoretical underpinning of this assumption. The physical nature of the gradient 
terms of our generalized micromagnetic model for ferromagnets with softened longitudinal 
fluctuations is explained. As part of the discussion we also address the relationship of our 
generalized micromagnetic model and the spin-fluctuation theory of weak itinerant-electron 
magnets. In section lTTD we review the experimental evidence for r\ < 1 in real materials. This 
is followed by a discussion in section |TV| of available evidence that supports a spontaneous 
skyrmion ground state in the itinerant-electron magnet MnSi. The supplement concludes 
with a discussion of the broader implications of spontaneous skyrmion ground states in 
magnetic metals in section IVl 



II. THEORETICAL UNDERPINNING OF THE LONGITUDINAL STIFFNESS 

In the following we address from a theoretical point of view the possiblty of reductions 
of the longitudinal stiffness in magnetic materials and their importance. For this purpose 
we consider the ratio rj of the longitudinal to the transverse magnetic stiffness, where 77 is 
identical to that in Eqn. (1) in the main text In our discussion we focus on the ordered 
state of isotropic ferromagnets. However, the arguments can be generalised in principle 
to all magnetic materials. Before entering a more technical level we present an intuitive 
explanation. 

It is instructive to consider 7] at first for isotropic local moment ferromagnets (LMFM). 
The conventional description of LMFM leads to a continuum theory for the magnetization 
distribution with 77 = 1, which is derived as the continuum limit of a Heisenberg model 
(for an elementary derivation see, e.g., Q). As shown by Patashinskii and Pokrovskii Q], 
magnon-magnon scattering in the ordered state forces the longitudinal stiffness to vanish 
for vanishing wavevector q. This effect originates alone in the broken continuous symmetry 
of the isotropic Heisenberg magnet, where rj = 1 for q — * (see also, Zwerger j^J). On more 
general grounds we note that the conservation of total spin, based on the SU(2) symmetry 
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of the spin system, requires that a reduction of magnetisation amplitude has to occur as 
a higher order process that involves at least two transverse processes. The finite lifetime 
of transverse excitations in turn results in a reduced longitudinal stiffness. If there are 
microscopic processes limiting the lifetime of transverse excitations in addition to magnon- 
magnon scattering, the longitudinal stiffness should be reduced further. Examples for such 
processes include: (i) interactions of magnons with either intraband or interband particle- 
hole excitations in itinerant magnets, (ii) interactions of magnons with defects, and (iii) the 
interaction of magnons with dipolar components of the magnetisation. 

The stiffness of a magnetic state is a consequence of the correlations introduced by these 
various microscopic processes. In an ordered or field-polarized state, the correlation lengths 
transverse and longitudinal to the local polarization direction, £j_ and £y, respectively, are 
generally different. These correlation lengths are a direct measure of the stiffness. Because 
the decay of longitudinal excitations requires higher-order transverse processes, we expect 
that £|| is always shorter than Thus, the longitudinal stiffness is expected to be lower 
than the transverse stiffness, i.e., in real materials we expect 77 < 1. This leads us to conclude 
that the longitudinal stiffness in real materials is always reduced below pure magnon-magnon 
scattering (rj = 1), whence 77 < 1. This conclusion is strongly supported by the experimental 
data reviewed in section IHIl which gives in the local moment ferromagnet EuS (t/euS = 
0.925), and the itinerant-electron ferromagnets Ni (r/ Ni = 0.65) and MnSi (r/MnSi — 0.4). 

We now enter the more technical part of the discussion. It is helpful to clarify at first the 
terminology we use. When we refer to 'micromagnetic models' we refer to a language that is 
in the tradition of the continuum theory of magnetism founded mainly by Landau's school in 
the 1940s. Micromagnetic models are used in the field of technical magnetism and magnetic 
domains for length scales down to the nanometer scale. Traditionally micromagnetic models 
consider only changes of the magnetization direction, while the amplitude is fixed by the 
saturation magnetization. Our generalised micromagnetic model for the magnetic stiffness 
is hence outside the tradition of micromagnetism. The field-theory underlying conventional 
micromagnetic models is the non-linear a model (NLcrM). In contrast, the description de- 
veloped for spin fluctuations in itinerant-electron magnets is referred to as spin fluctuation 
theory 0, Q . Spin fluctuation theory is based on Fermi liquid theory and the spontaneous 
spin polarisation of the metallic state. Traditionally, it considers materials where the mag- 
netization is free to change direction and amplitude. Both, micromagnetic models and spin 
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fluctuation theory can be represented in terms of a Ginzburg-Landau functional of a 3- vector 
model for the magnetization distribution. 

The task of this discussion is to bring out how to modify micromagnetic models for a 
description of itinerant-electron magnets. The discussion this way also aims to clarify the 
nature of the gradient term in 77. In the following we first consider the relationship of 
the conventional action of spin fluctuation theory with the NLcjM as the low temperature 
field-theoretical approximation that corresponds to micromagnetic models. This allows to 
illustrate how to change micromagnetic models to become appropriate for itinerant-electron 
magnets. We then provide arguments that mechanisms beyond magnon-magnon processes, 
which additionally reduce the longitudinal stiffness, indeed lead to a significantly reduced 
ratio 77 < 1 in real materials. 

We begin with an outline how to derive a corresponding micromagnetic model from spin- 
fluctuation theories for weak itinerant-electron ferromagnets. The appropriate Ginzburg- 
Landau functional in spin fluctuation theory is developed in a classical vector m(x) with 

n q 

three-components where the action (or reduced Hamiltonian) |2J is given by 

S(m) = j d d x[A(^2d a mf + r|m| 2 + 6|m| 4 ...]. (1) 

J OL 

The parameters b and r = a(T — Tc) are expansion coefficients with Tc the Curie tempera- 
ture. The parameter A is the phenomenological exchange stiffness parameter. 

To establish the relationship with the micromagnetic model we transform the action into 
a form reminiscent of the NLaM by writing m(x) = m(x)n(x) with |n(x)| = 1. Here 
m(x) describes the amplitude of magnetization, while n(x) is a unit vector that accounts for 
changes of the direction. In Equation (0) we now replace m(x) = m(x)n(x). This generates 
two gradient terms in the action 

5(771,11)= j d d x[Am 2 {x)(^d a nf + A^(d a m(x)) 2 + rm(x) 2 + bm i ...}. (2) 

OL OL 

Eqn. (J2J reduces to the NLaM when the amplitude is kept constant, m(x) = m = constant, 
i.e., the fluctuations of the modulus m are neglected. It follows 

S (n) = j d d x[Am 2 (J29 a n)% (3) 

OL 

where constant terms have been dropped. This is the low-temperature field-theoretical 
approximation used in micromagnetism. To go beyond this approximation, we have to 



include fluctuations in m(x). Considering that the longitudinal and transverse stiffness may 
differ in the ordered state, we use two different parameters for the two gradient terms in 
Eqn. (J2J). This is achieved by introducing the parameter 77. 

Si(m,n) = / d d x [Am 2 (x] (^<9 Q n) 2 + Ar\ ^(<9 a m(x)) 2 + rm(x) 2 + bm(x) 4 . . .] . (4) 

a a 

As shown in the following, rj measures the ratio of longitudinal to transverse stiffness. To 
bring out the nature of the gradient terms in the context of our micromagnetic Ansatz, 
Eqn. (J3J) can be rewritten as a functional of the magnetization vector m and amplitude m 

Si(m,m) = / d d x[A{J2d a m) 2 - {I - rj) A J](9 a m(i)) 2 + rm(i) 2 + 6m 4 ...]. (5) 

a a 

For rj < 1, the form of the action Si(m, m) implies that fluctuations of the amplitude lead 
to energy gains as compared to the case rj = 1. In other words, the action Si(m, n) (EqnHJ) 
describes softened longitudinal stiffness for rj < 1. 

The fixed modulus version of the action given in Eqn. which may be derived by 
approaching the low-temperature limit by taking r — > —00, b — ► 00, and r/b =const, has 
been analysed earlier in the context of critical phenomena and phase transitions 0, |t| . In 
the present context, we are interested in the gradient energy within the magnetically ordered 
state at an intermediate length, which is set by the competition between microscopic direct 
exchange and the chiral Dzyaloshinskii-Moriya exchange. It has been shown for the long- 
distance behaviour that only the usual gradient energy describing the rotation of a fixed- 
length magnetisation is relevant as used in micromagnetic theory in Eqn. © ^, 0, Q- 
More specifically, in the renormalization group analysis one finds that the gradient term in 
(1 — rj) in Eqn. © is irrelevant for the critical properties of ferromagnets i.e., in the 
limit of diverging correlation length at Tc- The analysis moreover stresses that only the 
transverse fluctuations are relevant to drive the phase-transition from the ordered state into 
the isotropic paramagnetic state. At the same time it is found that the irrelevant term with 
factor (1 — rj) has a different scaling behaviour as compared to the transverse stiffness. 

nn 

These renormalization group results [8|, |9( assert that the generalised gradient energy in 
our model does not change the expected universal behaviour of an isotropic ferromagnet 
with collinear ordering. This is the formal proof that we may freely use a generalized 
version of the gradient energy with rj < 1. As the two gradient terms have differing scaling 
behaviour, a coarse-grained effective action after a suitable renormalization procedure will 



always lead to gradient terms with 77 7^ 1 for the ratio of effective longitudinal and transverse 
stiffness. Therefore, this ratio will also depend on the scale considered under coarse-graining. 
When taken together, these theoretical results justify the asssumption that we can represent 
the two different longitudinal and transverse correlations lengths in the ordered state by 
the generalization of the gradient energy, which is embodied in the two parameters A and 
77. For arbitrary values of 77 > the gross qualitative properties of the paramagnetic to 
ferromagnetic transition and the ordered state are unchanged. However, the term arising for 
i] < 1 in the phenomenological theory given by Eqn. (J3J is the most important and simplest 
term to account for effects due to two different lengths for the stiffness of the ordered state. 
This term changes only non-universal properties. Essentially, it leads only to a shifted 
ordering temperature |8|. Therefore, the effects due to a softened longitudinal stiffness for 
rj < 1 are difficult to trace in the properties of a conventional collinear ferromagnet. However, 
the generalization for rj < 1 becomes important at intermediate lengths in the ordered state 
when competing magnetic couplings are present. In the case of chiral magnets considered 
in Ref. Q] a twisting length arises due to the competition between Dzyaloshinskii-Moriya 
couplings and direct exchange. This gives an additional length scale. The chiral couplings 
have a tendency to prompt the formation of inhomogeneous multiply twisted states. This 
is associated with the shift of the ordering temperature, according to the temperature scale 
T D as defined in Eqn. (3) of Ref. 

For the resulting multiply twisted non-collinear magnetization distribution, the difference 
of longitudinal and transverse stiffness becomes crucial. The formation of inhomogeneous 
non-collinear magnetization structures is faciliated by local shifts of the transition temper- 
ature that are enabled through a softened longitudinal stiffness. This is the reason why the 
generalized gradient term of our model with 77 < 1 has to be taken into account. We note 
that the term we consider is the leading-order correction needed to generalize the gradient 
term in a continuum model for the magnetization. Thus the conclusions of our study are 
intuitive and of the greatest generality possible. 

While our work was going on, it has been suggest ed that higher order gradient terms 



may also stabilise multiply twisted structures |11L Il2j. However, the possible existence of 
these terms is specific to certain materials and difficult to justify microscopically, so that 
the conclusions are not physically transparent and lack general relevance. 

Having justified the gradient term in 77 purely from a technical point of view, we now 
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consider how this term evolves under coarse graining in a general context. This serves to 
certify that the mechanisms listed above, which in principle may reduce the longitudinal 
stiffness, indeed will generate an effective model with t] < 1. Following convention, a gen- 
eralised Ginzburg-Landau theory may be written in terms of transverse and longitudinal 
fluctuating fields, 4>±{x) and m(x), respectively [ljj. The longitudinal component or modu- 
lus may be further decomposed as m(i) = M + m'(x) with a constant expectation value of 
the magnetization M and fluctuations of amplitude m'(x). This leads to the action 

S e (m',<f> ± ) = I d d x [A (^2d a (f>±) 2 + Ar]^2(d a m'(x)) 2 

a a 

+ rM 2 + bM 4 ... 

+ r\\ m'(x) 2 + rj_ |0±(x)| 2 

+wi m'(x) \<t>±(x) | 2 + w 2 m'(xf + + b (m' 2 (x) + \</>±(x) | 2 ) 4 ... (6) 

with constants r\\ = r + 32 bM 2 , r± = r + 8bM 2 w± = w 2 = 16 bM 2 . We note that the action 
given in Eqn. © is formally the analogue to that considered for itinerant-electron magnets 
at the border of ferromagnetic order p. The limit of a self-consistent coupling of finite q 
modes to the uniform magnetisation (q = 0) of this model has been studied extensively by 

n n 

Moriya |6J, Lonzarich |5[ and others. An important assumption of their analysis is that the 
response of the system is that of a Fermi liquid in the random phase approximation (RPA). 

To establish the relationship of our micromagnetic model, given by Eqn. (1) in the 
manuscript 1] with Eqn. (jUJ) it is in principle necessary to carry out a renormalisation 
of the action S e (m', 0jj up to the length scales addressed in the micromagnetic model. Un- 
fortunately, a rigorous renormalization analysis below the ordering temperature for isotropic 
n-vector models has not been achieved due to great technical difficulties in treating mode- 
,„ode coupling co„,c, y Q Q. TKe 8 e dirties may be t raCed to the eritieality of the 
transverse field components along the whole coexistence line below Tc and the associated 
system of Goldstone- modes 0, 13 16 , 17 1 . The apparent non-analytic form of the gradient 
energy for rj < 1, as seen in Eqn. (jl), is likewise expected because of the divergencies related 
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12|. 



to transverse fluctuations in the ordered state 

While it appears a major challenge to perform a rigorous renormalisation of isotropic n- 
vector models in the ordered state that is well beyond the scope of the work reported here, 
it is nevertheless possible to track the evolution of rj under renormalisation in a Gedanken 
calculation. We begin by noting that in a perturbation expansion beyond one-loop-order, 



the renormalization of the fields <j)± and m' scales as powers of the correlation lengths 
and £||" d<1> , respectively, where the exponent d$ is also known as anomalous dimension of the 
fields Il8l. This implies that renormalizations would, in principle, break the isotropy of the 
model 13[- To guarantee the isotropy of the model under renormalisation, the longitudinal 
and transverse field components should instead be rescaled consistently by the same factor 
Thus the action has to be re-paramaterised under renormalisation, which is achieved by 
forcing the term in rj to receive a renormalization by an additional factor {C\\/C±) di - We note, 
that along or near the co-existence line, the anomalous dimension c?$ is to be considered 
as an effective exponent of a cross-over region Q]. Therefore, it will depend on distance 
to the critical points at temperatures T = Tc and at T = 0. Because the longitudinal 
correlation length on general grounds is shorter than the transverse correlation length, as 
discussed above, rj must remain less than one also after integrating out degrees of freedom 
for the construction of effective coarse-grained theories. 

A rigorous proof along the lines of our argument is difficult, because the transverse cor- 
relation length diverges in zero magnetic field, £j_ — > oo, due to the system of transverse 
Goldstone modes. To avoid the problems stemming from the divergence of the correlation 
lengths, we may assume that the divergencies at long lengths are cut-off by inclusion of an 
arbitrary small field in longitudinal direction. This field serves also to define the longitu- 
dinal and transverse directions. Then, the scaling properties under coarse-graining can be 
connected with spin-fluctuation theories of an isotropic itinerant electron ferromagnet in the 
presence of finite polarization, which generally also yields different perpendicular and lon- 



gitudinal correlation le ngth s as discussed in Ref. [19|, 20]. The form of the spin-fluctuation 
theory reported in Ref. [3, Q takes mode-mode coupling into account, and the approach 
assumes that the squared modulus of spin-fluctuations is held constant. Thus, it can be 
viewed as an approximation for an itinerant-electron ferromagnet with maximal longitudi- 
nal rigidness, similar to the case of a Heisenberg ferromagnet with local moments (constant 



modulus). From Refs. 19|, |20j, |2l| we have that < 1 everywhere below the ordering 



temperature at finite polarization, and — 1) asymptotically at Tq- The relation 

(£||/£±) < 1 given from the spin- fluctuation theory of an itinerant-electron ferromagnet in 
the ordered state is the expected result from our general consideration on transverse and 
longitudinal correlation lengths. Consequently, following our scaling argument, 77 < 1 must 
hold below the transition temperature as a leading order (irrelevant) correction in a coarse- 
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grained theory for the itinerant-electron ferromagnets. We note, that this result does not 
depend on the detailed electronic structure, say, whether a single-band or multiband system 
is considered. 

In summary we conclude that the assumption rj < 1 should hold quite generally in effective 
micromagnetic theories, for the ordered magnetic state whenever amplitude fluctuations of 
the magnetisation are important. For fixed moment magnets that are well described on a 
microscopic scale by a Heisenberg model, the effect may be rather weak as it involves higher- 
order perturbation contributions. If the underlying microscopic processes in a material 
allow for longitudinal fluctuations on short lengths scales, an effective term 77 < 1 is always 
expected. 



III. THE LONGITUDINAL STIFFNESS IN REAL MATERIALS 



In the following we review the transverse and longitudinal magnetic stiffness observed 
experimentally by polarised neutron scattering at the border of magnetic order in real ma- 
terials. In the context of our micromagnetic model we are interested, if the longitudinal 
stiffness is reduced on intermediate lengthscales. For this purpose we consider the ratio rj 
of the longitudinal to the transverse magnetic stiffness, where we refer to section |H] for a 
more detailed theoretical underpinning of the analysis. As materials we consider the cubic 
systems EuS, Nickel and MnSi, which exemplify local-moment ferromagnetism, itinerant- 
electron ferromagnetism and weak itinerant-electron ferromagnetism, respectively. In MnSi 
we the presence of weak chiral interactions do not alter the outcome of these estimates. 

In the three compounds considered here the wavevector dependent susceptibility follows 
the Ornstein-Zernicke form Xi l {q) ~ c i{ K 1 + q 2 ) + ■ ■ ■ over major portions of the Brillouin 
zone, where the index i denotes transverse (i =_L) and longitudinal components (i = ||), 
respectively. This form is consistent with the generalised spin fluctuation model given in 
Eqn. ©. The parameter Ki = is the inverse correlation length and the parameter Cj the 
stiffness. The stiffness ratio 5 = cu/c± permits a comparison of different materials without 
need for complicated normalisation procedures of the data. For the cubic systems considered 
here the total susceptbility near Tc is expected to be consistent with a 3-component vector 



model 



22j . where Xtot — 2x_i_ + Xll accounts for 2 transverse and 1 longitudinal mode. For a 



total susceptibility of the Ornstein-Zernicke form it follows that Xtot = c tot(^t t + 9 2 ); where 
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transverse and longitudinal contributions are given by x± — C ±( K ± + <? 2 ) = {3/2)Xtot an d 
X^ 1 = C||(/Cm + g 2 ) = 3xtot, respectively. Thus, when the longitudinal and transverse stiffness 
are the same it follows 5 = 2. The parameter 77 introduced above is given by rj = 5/2, i.e., 
77 is the stiffness ratio per mode. 

By definition x(o) is equivalent to the frequency integrated scattering intensity. To obtain 
the stiffness ratio of EuS we have taken the slo pes of the inverse integrated scattering 
intensities, and thus xT 1 ^)) shown in Fig. 6 in Ref . |23l| . Here we obtain <5e u s — c ||/ c ± ~ 1.85. 
This corresponds to t^euS ~ 0.925. Taking likewise the ratio of the slopes of the integrated 
scattering intensity for the itinerant-electron ferromagnet Nickel shown in Fig. 10 in Ref. [2^ 
we obtain 5^i ~ 1.3 or 77™ ~ 0.65 Finally, for MnSi the available spin-polarised neutron data 
yield ^MnSi ~ 0.8 24j,|25j], i.e., ?7MnSi ~ 0.4. Hence, in the range where data are available, 
the longitudinal stiffness in Ni and MnSi is reduced by 35% and 60%, respectively, Because 
the g-range considered in Ni and MnSi is in a regime where neutron scattering shows an 
abundance of spin-flip particle-hole excitations, we conclude that particle-hole excitations 
are dominant in causing the reduction of the longitudinal stiffness of Ni and MnSi. 

The conjecture that particle-hole excitations are dominant in reducing 77 is supported 
further by the characteristic frequency dependence, which reflects on the mechanisms that 
determine the transverse and longitudinal stiffness. For local moment systems such as EuS, 
it is given by the Resibois-Piettes function. In contrast for itinerant-electron magnets, an 
additional softening has been noticed (see, e.g., Fig. 13 in Ref. Q for Ni). A simple 
measure of the frequency dependence and thus 77 may be the spin fluctuation temperature, 
T s f, defined from the width of the spin fluctuation relaxation frequency spectrum. The 
ratio of Curie temperature to spin fluctuation temperature, r = Tc/T s f, may be used as 
an estimate for the magnitude of 77, notably, the smaller the ratio r the smalller 77. This 
suggestion is consistent with the experimentally observed reduced longitudinal stiffness in 
Ni and MnSi, discussed above. 

When taken together, the experimental evidence establishes beyond doubt that many 
materials exist, which may be described by the regime 77 < 1 of our model. Moreover, the 
experimental data identifies MnSi, which supports chiral interactions, as a prime candidate 
for a spontaneous skyrmion ground state, as discussed in section HVl Interestingly the exper- 
imental data even shows 77 < 1 in the local moment magnet EuS. This is not a suprise, since 
a multitude of mechanism exists that reduce the longitudinal stiffness. As a consequence, 
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spontaneous skyrmion phases may exist in chiral local moment magnets, albeit they may be 
much less stable. 



IV. POSSIBLE EVIDENCE OF SYKRMION GROUND STATES IN MNSI 

As the main result of our study we predict the spontaneous formation of a skyrmion 
ground state in a temperature interval prior to helical order in all materials with soft mag- 
netisation amplitude and weak chiral interactions. Bulk compounds which develop long- 



wavelength helica 



i Q, 



FeGe 



order of the kind considered here are, for instance, MnSi 
and (Fe,Co)Si [29 1, where the modulation wavelengths are AMnSi ~ 180A, AfcGc ~ 700A 
and A(F e ,co)Si ~ 300A, respectively. In these materials the helical modulation is driven by 
Dzyaloshinsky-Moriya chiral interactions that originate in the lack of inversion symmetry in 
their cubic B20 crystal structure. Being transition metal compounds the magnetic order is 
due to itinerant electrons, although this has ou.y been shown experiments in MnSi Q. 
As discussed above polarised neutron scattering data exist for MnSi [24J that establish a 
reduced longitudinal stiffness (rj MnSi = 0.4). Thus at least MnSi meets all the requirements 
for the formation of the predicted spontaneous skyrmion ground states. 

The helical order in MnSi, FeGe and (Fe,Co)Si may be seen most easily in small-angle 
neutron scattering (SANS) experiments, where sharp magnetic satellites were observed. In 
a small temperature interval above the onset of helical order, SANS data for MnSi and FeGe 
show neutron scattering intensity located on the surface of a sphere in reciprocal space as 
an additional feature^. I3I | - 

The radius of this sphere corresponds to the same modulation length as the helical order. 
The sphere of scattering intensity has therefore traditionally been interpreted as fluctuating 
helical order. However, in these studies it was not considered that the scattering intensity 
may also be related to a frozen-in form of order. In this respect it is interesting to note, 
that the scattering intensity on the surface of the sphere is perfectly consistent with an 
amorphous array of randomly oriented cylindrical skyrmion tubes. 

To better identify of the ring of scattering intensity we present a rough estimate of the 
temperature interval AT = To — Tc, for which the skyrmion ground state is expected in 
MnSi. The lattice constant of MnSi is aMnSi — 

4.55 A. There are four formula units per 
unit cell, so that the number-density is N/V = n = 4.246 ■ 10 28 m~ 3 . The effective moment 
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determined from the Curie- Weiss temperature dependence of the susceptibility is p e g = 2.2/j,b 
(3^] and the susceptibility \ — Ccw/{T — Tc) where C*cw = nfiofi B (jp%ff) / (3kB) = 0.534 K. 
The upper limit of the skyrmion transition temperature is T^ = Tq + (1/2)(D / (2A))(D / a) 
as given in Eq.(2) in Ref. Q]. These parameters have to be taken from experimental data, 
notably Ref. [33j . The parameter a is given by the initial susceptibility and may be estimated 
from the Curie- Weiss susceptibility: a = l/(2Ccw) — 0.936 K -1 . The modulus of the helical 
propagation vector is go = (D/(2A)) = 0.039 A" 1 . In the following we use temperature 
as unit of energy- density. The exchange stiffness A may then be estimated from A ~ 
(c4 nSi )T c = 50 • 10" 3 eV A 2 = 587 K A 2 . The strength of the DM-coupling D may be 
estimated from q as D = 1.9 • 10 -3 eV A= 22 K A. These parameters finally result in 
T h = T c + (l/2)q (D/a) = T c + q DC C w = (29.5 + 0.46) K. The ordering temperature of 
the skyrmion phase (according to Eq.(3) in Ref. Q]) is obtained by doubling the difference 
between Tc and so that AT = Tp — Tc = 0.9K. This value of AT gives an order-of- 
magnitude estimate for the temperature interval of the stable skyrmion phase in MnSi as 
shown in the phase diagram of Fig. 2 in our manuscript bj . 

Detailed studies of the specific heat (Fig. 6 in Ref.j^J, 3]) and thermal expansion [2^ 
of the paramagnetic to helical transition in MnSi show the presence of a spike at T^ and 
a broad shoulder in a temperature interval of order ~1 K above this transition. Detailed 
SANS studies of the ring of intensity as function of temperature suggest that the ring 
extends over the same temperature interval where the shoulder in the specific heat is seen 
Interestingly the temperature interval corresponds fairly well with that predicted for 
the skyrmion phase. In combination, the broad shoulder in the specific heat clearly signals 
the presence of quasi-static order already above T h that may be either related to randomly 
oriented helical domains or an amor pho us texture of skyrmions. An important additional 
aspect that has been emphasized in 34j is, that the shape of the specific heat anomaly in 
very high purity samples is quite different, yet, still supports the possible existence of a 
precursor phase. However the shoulder is observed for samples with largely varying residual 
resistivity ratio |25j. It is possible that the dependence on extreme sample purity signals 
that the stabilisation of the quasistatic order sensitively depends on residual defects in the 
material. 

The spontaneous formation of skyrmion ground states allows to explain many of the 
mysterious properties of MnSi at high pressure and low temperatures. It was recently 
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reported by one of us (CP and collaborators), that the metallic state of MnSi changes as 
function of pressure from a quadratic temperature dependence of a weakly spin-polarized 
ermi liquid to a T 3 / 2 resisitivity that is stable over a remarkably large pressure range 



38, 



39 



spects appear to be inconsistent with the present day understanding 
of metallic magnetism: First, as pointed out in Ref. [38j a T 3 / 2 resistivity is normally observed 
in spin glasses and amorphous ferromagnets, where it is explained by a diffusive motion of 
the charge carriers 42]. However, for the high-purity MnSi single-crystals investigated the 
behaviour of a moderately enhanced Fermi liquid is expected by all accounts. It was therefore 
emphazised that the diffusive charge carrier motion would have to be intrinsic. Second, the 
T 3 / 2 resistivity is extremely stable as function of pressure, where it has been seen up to at 
least 45kbar^ 3p c . This suggests the formation of an extended new phase in contrast to 
a cross-over phenomenon normally observed at quantum phase transitions. Third, the T 3//2 
resistivity is extremely stable as function of magnetic field, where experiment shows that it 
collapses abruptly only above a certain critical field related to a so-called itinerant-electron 



metamagnetic transition (MMT) |43j. At the highest pressures the T 3 / 2 resistivity extends 
at least up to IT (cf Fig. 3 in Ref. (^J) which substantially exceeds the field needed to 
collapse helical order B c « 0.6T. We note that B c studied by various magnetic probes is 
insensitive to pressure. 

As a possible resolution to these contradictions we suggest that the Fermi liquid to non- 
Fermi liquid transition in MnSi may signal a transition between helical order and an amor- 
phous skyrmion ground state. A spontaneous amorphous skyrmion ground state provides 
a simple intrinsic mechanism that explains why the transport properties of an amorphous 
ferromagnet are seen even in a high purity material in the absence of disorder or frustration. 
The leading order effect of pressure in this scenario would be to generate a slight mode soft- 
ening that reduces the stability of helical order, but not to collapse the magnetic moment 
on local scales completely as normally assumed in quantum phase transitions. Because the 
skyrmions form as a genuine, stable ground state, it would finally not be difficult to explain, 
a priori, the extent of the T 3//2 resistivity behaviour in the temperature-pressure phase- 
diagram. Also, it has been shown that the core of cylindrical skyrmion tubes is extremely 
stable against external fields Q|. Thus, if the T 3 / 2 resistivity is due to an amorphous 
skyrmion texture it is expected that the behaviour survies to much higher fields. The signa- 
ture of such an amorphous skyrmion ground state in neutron scattering would be intensity 
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on the surface of a small sphere in reciprocal space, similar to that seen at ambient pressure 
above Th- 

Recent neutron scattering studies of the magnetism as function of pressure indeed show 
that large ordered moments survive deep into the phase where the T 3 / 2 resistivity is observed 
The ordered moments are organised such that they lead to scattering intensity on the 
surface of a sphere in reciprocal space. There is, however, a very important difference with 
the intensity on a sphere seen at ambient pressure near Th. The ambient pressure intensity is 
isotropic, while broad maxima are observed for (110) at high pressure. The (110) direction 
is unusual for this cubic magnet. This is so, because (i) the leading magnetocrystalline 
anisotropy favours easy magnetic axes only in (111) or (100) directions. And, (ii), in itinerant 
d-electron magnets as MnSi with very weak anisotropy, it is very unlikely that higher-order 
anisotropy contributions do play any role. In fact, the easy axis of the helical order in MnSi is 
(111) and in FeGe (100). Thus, if there was randomly oriented helical order at high pressure 
in MnSi, it would be very difficult to explain the broad maxima along (110). In contrast, for 
a condensate of skyrmion tubes the odd maxima for (110) are easily explained, if we suppose 
that the easy magnetic axis of MnSi, (111), does not change as function of pressure, while 
the magnetism changes from a helical modulation along (111) to an amorphous texture 
with cylindrical skyrmion tubes having their axes preferentially along (111). For cylindrical 
skyrmion tubes that are stratified along (111) neutron scattering intensity in reciprocal space 
is expected on great circles perpendicular to (111) directions. As all (111) directions produce 
intensities along such great circles, these will intersect and produce maxima of intensity in 
the (110) directions. The available experimental data are consistent with this possibility. 

However, in contrast to the remarkable extent of the T 3//2 resistivity under pressure and 
magnetic field, the sphere of scattering intensity under pressure is seen only in a fairly 
small portion of the phase diagram. A possible explanation discussed in Ref. |45j is that 
the textures are dynamically destabilized and become 'invisible' in the neutron scattering 
experiment. Such a slowly meandering form of magnetic order could persist on a time scale 
that it still much slower than the time scale relevant to the electrical resistivity. A liquid 
instead of a frozen skyrmion phase in this pressure range provides a natural explanation for 
the difference in the transport measurements and the visible static magnetic order. Such 
a liquid skyrmion phase would not be unusual and may be analoguous to liquid vortex 
phases in superconductors |4a|. It will, in any case, require a host of novel experimental 
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techniques to settle this issue unambiguously. For instance, real space imaging of helical 
order in thin (Fe,Co)Si layers was recently reported using Lorentz force microscopy (LFM) 
|47|. However, samples that can be studied with LFM have to be atomically thin so that 
uncontrolled mechanical strains and their interplay with the magnetism represents a major 
concern. As regards LFM in MnSi we note, moreover, that even state-of-the-art LFM does 
not offer sufficient resolution yet to image the magnetic structure with sizes of the twisting 
length 18 nm. We believe that, as an important aternative to real-space imaging, the most 
promising method in the immediate future will be full three-dimensional polarisation analysis 
for neutron scattering. This is presently under development under the synonym mu-PAD. 



V. ON THE BROADER IMPLICATIONS OF SPONTANEOUS SKYRMION 
GROUND STATES DUE TO CHIRAL INTERACTIONS 



Skyrmion states in field theories have attracted great interest because they represent 
stable, particle-like objects with a nontrivial topology that originates in fundamental sym- 



metries (48J, 49]. For instance, Skyrme's field theory for elementary nuclear particles has 
been applied to condensates of a few skyrmions and extended lattice structures 51] as 
models of ordinary nuclei and compact stars, respectively. Our prediction of spontaneous 
skyrmion ground states in chiral magnets is connected to the fundamental questions about 
the relation between continuum theories, the countable nature of particles and the possibility 
to relate both aspects by topology in extended matter systems. In this more general context 
he role of chirality and chiral interactions, which are possible in condensed matter systems 
5^ | , has not yet been appreciated. Chiral interactions in condensed matter systems offer an 
exceptionally rich setting, because they exist in many different contexts, e.g., (i) spin-orbit 
interactions in non-centrosymmetric magnetic materials, also referred to as Dzyaloshinsky- 

54], (ii) in non-centrosymmetric ferroelectrics J55, 0], (iii) 



Moriya (DM) interactions 



53, 



for certain structural phase transitions 



57 



581 ] . (iv) in chiral liquid crystals 
60,0. 



or (v) in 



the form of Chern-Simons terms in gauge field theories 

The example of an extended magnetic skyrmion texture shown in Figure 4 of Ref. Q] 
demonstrates the condensation of particle- or string-like localized objects. There are, 
however, further intriguing possibilities beyond these extended ground-states. Magnetic 
skyrmions may also arise as non-linear localized excitations both in homogeneously ordered 
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systems |44] and in the paramagnetic state. This possibility has not yet been considered in 
theories on the metallic state of systems with broken inversion symmetry. In particular, the 
description of skyrmion excitations on a paramagnetic background would be very similar 
to that of particles generated in a vacuum state. While being large and localised, they are 
subject to an intrinsic length scale related to the twisting length of the system. This length- 
scale implies non-dispersive behaviour unlike propagating waves. Moreover, skyrmions in 
condensed matter systems may allow to study processes like particle generation and annihi- 
lation, scattering between particles, or particle confinement. 

Fundamental processes of this kind are usually the subject of field-theories in particle 
physics only. Images of such processes for elementary particles can therefore be inferred at 
best only indirectly. In chiral magnetic condensed matter, these fundamental processes take 
place in a continuum distribution of magnetization that may be visualized in a solid state 
laboratory by modern magnetic real- and reciprocal space imaging methods. In fact, because 
the spontaneous skyrmion ground states predicted here represent large (several hundreds A) , 
but localized non-linear spin-fluctuations, they would provide the first cases of skyrmions 
in a condensed matter system, where detailed imaging of such processes may be possible. 
With this perspective in mind, it is gratifying that novel chiral crystals have recently been 
found, where superconductivity coexists with magnetic order, as in the novel heavy-fermion 
superconductors CePt 3 Si 6^] and Ulr 0|. In these materials different types of long-range 
order coexist, so that it may be possible to study interactions between skyrmionic excitations 
in the associated order-parameter fields which in turn have different symmetries. In fact, 
chiral condensed matter systems could eventually be used to mimick a universe with different 
types of particles and their interactions. 

In spite of their low symmetry, metallic magnets with broken inversion symmetry are 
rather common. In particular, Dzvaloshinsky-Moriya interactions are induced by surfaces 
in all magnetic nanostructures Thus magnetic nanostructures can endow science 

with an arena, where creation, interactions, and condensation of particle-like states in a 
continuous field can be engineered for experiment. 
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